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Abstract 

For a wide class of Hamiltonians, a novel method to obtain lower 
and upper bounds for the lowest energy is presented. Unlike perturba- 
tive or variational techniques, this method does not involve the com- 
putation of any integral (a normalisation factor or a matrix element). 
It just requires the determination of the absolute minimum and max- 
imum in the whole configuration space of the local energy associated 
with a normalisable trial function (the calculation of the norm is not 
needed). After a general introduction, the method is applied to three 
non-integrable systems: the asymmetric annular billiard, the many- 
body spinless Coulombian problem, the hydrogen atom in a constant 
and uniform magnetic field. Being more sensitive than the variational 
methods to any local perturbation of the trial function, this method 
can used to systematically improve the energy bounds with a local 
skilled analysis; an algorithm relying on this method can therefore be 
constructed and an explicit example for a one-dimensional problem is 
given. 

PACS : 03.65.Db, 05.45.Mt, 02.30.Tb, 02.60.Gf. 



1 Introduction 

In a large variety of interesting physical problems, finding the discrete spec- 
trum of an operator can be done with approximate methods only. Moreover, 
in most cases, it is a rather difficult task to estimate errors. Perturbative 
techniques often lead to non-convergent series and an evaluation of the dis- 
crepancies with the exact result is usually beyond their scope. For semi- 
bounded operators, variational methods naturally provide upper bounds for 
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the lowest eigenvalue and require much more work for providing a lower 
bound with Temple-like methods (Reed & Simon, 1978, XIII. 2). Another 
source of difficulties when dealing with perturbative and/or variational tech- 
niques, is that they both involve the calculation of integrals on the configu- 
ration space Q: the norm of the wavefunctions and some matrix elements of 
operators. This reflects the fact that the discrete spectrum of a differential 
operator encapsulates some global information on Q within the boundary 
conditions imposed on the normalisable wavefunctions. In this article, I want 
to propose an approximate method that will overcome these two obstacles: 
it can rigorously provide both lower and upper bounds without any kind of 
integration. Like the variational techniques, it will involve a (set of) trial 
normalisable function(s) with the appropriate boundary conditions and will 
concern in practice the lowest eigenvalue only. The bounds are given by the 
absolute extrema of a function defined on Q (the so-called local energy) . In 
a sense, this method allows to stay as local as possible in the configuration 
space Q: in order to improve the bounds, a local analysis near the extrema 
(or near the possible singular points) of the local energy is sufficient and 
necessary. This paper is organised as follows. In section |21 I give the proof 
of the inequalities that will be the starting point of the differential method. 
A comparison with what already exists in the litterature follows and the 
guidelines of the method are presented. Sections El lU and show how the 
method can be applied to three non-integrable quantum systems. Before the 
concluding remarks, in sectional I explicitely show on the quartic oscillator 
how the sensitivity of the differential method to local perturbations of the 
trial functions can be exploited to systematically improve the bounds of the 
ground state energy with an elementary algorithm. 

2 Bounding the ground state energy with the local 
energy 

Let us start with a quantum system whose Hamiltonian K acts on the 
Hilbert space of functions defined on a configuration space Q. Let us sup- 
pose that II has an eigenstate |<I>o) associated with an element Eq of the 
discrete spectrum. For any state the hermiticity of H implies the iden- 
tity ( I (-f^ — -E'o )!'/') = 0- If we choose |(/?) such that its configuration space 
representation (/?(g) is a smooth real normalisable wavefunction, we obtain: 



The crucial positivity hypothesis is to assume that we can choose one eigen- 
state such that its eigenfunction <I>o remains real and positive or zero in the 
whole Q. This generally applies to the ground state for which it has been 
shown in many cases that it is strictly positive in the interior of Q (Reed & 
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Simon, 1978, XIII. 12). Then, the real part of (dJ involves an integrand that 
is a smooth, real function constructed from the real part Hr of the differen- 
tial operator H. Then there exists a g in Q such that ^Q{q){Hfi — EQ)(p{q) 
changes its sign. Therefore, 

3qe Q such that {Hr - Eo)(p{q) = . (2) 

Let us now introduce a function on Q that is known as the local energy ^ 



« ^ . (3) 



From condition @, we immediately obtain that for all smooth real and 
normalisable state ip, 

inf(^L1(cz))^ Eo ^ sup (4) 

Q Q 

Surprisingly, these two inequalities are sparsely known in the literature 
(Barnsley, 1978; Baumgartner, 1979; Thirring, 1979; Crandall & 
Reno, 1982; Schmutz, 1985) and always under some more restricted con- 
ditions (the upper bound is often missing). 

The original proof presented here links the inequalities to the non-nega- 
tivity of the ground state without referring to the detailed structure of the 
Hamiltonian. In particular, it does not require for the Hamiltonian to have 
the purely quadratic form 

where o is a definite positive real matrix and V a well-behaved potential. 
Inequalities Q still apply (with the appropriate definition © of the local 
energy) in the presence of a singular potential, when there is a magnetic field 
and for an infinite number of freedoms (like in the non-relativistic quantum 
field describing a BEC condensate). Besides, it is not required for ip to be 
nonvanishing. It simply says that where ip vanishes faster than Hji ip, one 
or both of the bounds can be infinite and therefore useless. 

The first form of inequalities (jU (with its two bounds) is due to Barta 
(Barta, 1937) and was derived for the fundamental vibration mode of an 
elastic membrane. Though Barta writes that his method will be generalised 
in subsequent publications, I was not able to find any extensions of his orig- 
inal work before an article of Duffin (Duffin, 1947) where a Schrodinger 



""^The usual motivation for introducing the local energy is just to roughly estimate the 
dispersion in energy obtained for an approximated eigenfunction. For instance, when using 
Monte-Carlo methods for computing expectation values. The derivation of inequality I0J 
shows that for the ground state this qualitative approach can be made rigourous. 
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operator of the form H = —A + V is considered. Duffin shows that the 
Dirichlet boundary conditions imposed on the trial function (p can be re- 
laxed but he loses the upper bound. One obtains the equalities in (jlj for 
a flat local energy i.e. for (/? = $oj hence we will try to work with a ip 
that mimics the exact ground state best. Therefore, generalising the Barta 
inequalities by increasing the size of the functional space of (/?'s can be irrele- 
vant. One should instead keep working with a restricted set of trial functions 
that respects some a priori known properties of $0i such as its boundary 
conditions, its symmetries and its positivity^. 

More precisely, we will explain in the last part of this paper that, once a 
if that bounds the local energy is found, it is expected that there is only a 
finite number of independent directions in the functional space along which 
the bounds can be improved. In the following we will actually deal with 
a finite dimensional submanifold of trial functions <px where A stands for a 
small number of control parameters varying in a control space C. Accord- 
ingly, the strategy is clear: for, say, obtaining an optimized lower bound 
we will try to find sup^^g^ (infggg (£'ioc(A; (7))) where Eioc{X;q) stands for 

El^^\q). As long as the extremal values of q ^ £"100 ( A; g) can be followed 
smoothly with A (in particular the Morse points are generically stable), the 
problem is reduced to local differential calculations in C x Q in the neigh- 
borhood of the critical points: adding to a trial function an infinitesimal 
perturbation that is localized far away from the extremal point does not af- 
fect the energy bounds. One recovers the global sensitivity of the eigenvalue 
problem because the critical points of the local energy generically bifurcate 
for finite variations of A (Poston & Stewart, 1978; Demazure, 2000) and 
can jump to other distant points when a degeneracy occurs (Arnold, 1984, 
chap. 10, especially fig. 50). 

In mathematical physics literature, Barta's inequalities are always con- 
sidered within the context of the billiards systems (Laplacian spectra on a 
Riemannian manifold), even in the most recent papers (for instance (Bessa 
& Montenegro, 2004)). As far as I could search, the most advanced ex- 
tension to other physical problems has been made (tentatively) by Barnsley 
(Barnsley, 1978) but, for the same reasons as Duffin's (Duffin, 1947), he 
systematically loses the upper bound. Besides, he acknowledges he is unable 
to produce any non-trivial bound for the Helium atom. 



technicality should be mentioned here: If one chooses the trial states \(p) such that 
for all q £ Q, (qlf) > except, perhaps, for the a priori known zeros of (g|<l>o), we can 
deal with systems where the configuration variable q includes some discrete parameter like 
a spin index (the somewhat loosely notation Hinp{q) must be understood as the real part 
of {q\H\ip) and dq is the measure on Q possibly having a continuous and/or a discrete 
part). Indeed, under the positivity hypothesis, from Q we deduce that there must be a 
couple {q,q') in such that <4>{q) > 0, ip{q') > 0, {HR-Eo)ip{q) > 0, (HR-EoMq) ^ 0. 
Therefore inequalities I^J remain valid. For instance, if Hr is a (possibly finite) matrix, 
the local energy consists of a discrete (finite) set of real numbers. 
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One can easily understand Barnsley's failure: With variational meth- 
ods, a very rough estimation of the exact ground state wavefunction can 
lead to a reasonably good agreement for Eq while a simple local perturba- 
tion of the exact wave function can even make the local energy unbounded. 
Therefore, at first sight, one could see the sensitivity of the local energy as 
a major drawback of the method: variational methods are more robust to 
local perturbations of the trial function. But this argument can be reversed : 
compared to the rigidity of the variational methods, the differential method 
offers the possibility to improve the estimations at low cost provided we are 
able to implement a skilled strategy (eliminating the singularities, control- 
ling the behavior at infinity with JWKB techniques, increasing the absolute 
minimums, etc.). In the following, I explicitly show in many non-trivial 
cases, that once we have this strategy in mind, we can obtain interesting re- 
sults for complex systems. For instance, not only we can improve Barnsley's 
trivial bound for the ground state energy of the Helium atom, but it will 
be shown in section 0] how this result generalizes to any number of Coulom- 
bian particles. The calculations can be made analytically with a surprising 
simplicity. 

As far as the upper bound is concerned, the variational method leads a 
priori to a better approximation than the differential method since, for any 
normalised function (p, 

[ ^*{q)H^{q) dq= [ Mq)\'El%) dq ^ sup {e\^M ■ (6) 

Nevertheless, being free of any integration, the absolute maximum of the 
local energy is a quantity that is more easily accessible to analytical or 
numerical computations than the average value of H. 

3 Application to billiards; the example of the 2d- 
annular billiard 

As a first illustration of the differential method, let us consider the problem 
of finding the lowest eigenvalue ol H = —A/2 in a connected finite region Q 
with the Dirichlet boundary conditions imposed on B = dQ. Suppose that 
the boundary B is given by an implicit smooth scalar equation of the form 
b{q) = while the interior of Q is defined to be the set of points q such that 
b(q) < 0. Then, trial functions can be taken of the form if = fb for any 
arbitrary smooth function / that does not vanish inside Q. The only possible 
singular points of g i— > -E'l^j(g) are located on B and can be removed if / is 
chosen with an appropriate behaviour in the neighborhood of B. Imposing 
this behaviour for / is a priori a simpler task than solving the eigenvalue 
problem on the global Q: one dimension has been spared since we have to 
deal with some local properties of / near B. For instance, by generalizing 
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Barta's trick (Barta, 1937), one can easily check by a simple equation 
counting that when Q is algebraic i.e. when 5 is a polynomial, provided that 
we choose / to be a polynomial of sufficiently high degree n whose zeros are 
outside Q, we can find a polynomial g of degree n — 2 such that A(/6) = gb. 
Therefore E^^^ = —g/{2f) is bounded and finite upper and lower bounds 
of £"0 can be found. Let us apply this method to the asymmetric annular 
billiard that is an elegant paradigmatic model in quantum chaos (BOHIGAS, 
Boose, Egydio de Carvalho & Marvulle, 1993). B is made of two 
circles of radius 1 and r < 1 whose centers are distant hy S < 1 — r. Q is 
the 2d-domain in between the circles. The simplest choice of trial function 
is to take if{x,y) = b{x, y) = [x^ + — r^] [(x — (5)^ + — 1] . One can check 
analytically that the lower bound of e\^^ is 

■ f (Z^\ = / 8[(x-V2)^ + y^-(l+r^)/4] 

For r = 3/4 and 6 = 0.1, a simple numerical computation shows that ((7|) 
is finite and obtained at (x, y) ~ (0.86, 0) and leads to Eq ^ 28.390 to be 
compared with the exact result £"0 — 42.94. As one could have expected 
with the rough trial function chosen above, the estimation is not very precise 
but the calculations required here to get this result are much simpler than 
the ones involved in a variational method (that provides the complementary 
upper bound 55.32 with the same test function) or by the exact numerical 
resolution that requires to find the smallest root of an infinite determinant 
made of Bessel functions. 




4 The many-body Coulombian problem 

The next examples, presented in this section and in the following, will illus- 
trate that the first strategy for obtaining finite bounds is to get rid of the 
singularities that may appear in the local energy. When Q is not bounded, 
one must have a control over the behaviour of the trial functions as q goes to 
infinity. For a multidimensional, non separable, Schrodinger Hamiltonian, a 
JWKB-like asymptotic expression is generally not available (Maslov & Fe- 
DORIUK, 1981, Introduction). Nevertheless, the differential method is less 
demanding than the semiclassical approximations: we will ask that the local 
energy be bounded at infinity but we will not require it to tend to the same 
limit in all directions. As already shown in the annular billiard problem, for 
the sake of simplicity one could start with a less ambitious program and try 
to obtain just one nontrivial inequality in ((1]). 

The second example is to consider a system of N non-relativistic, spin- 
less, charged particles living in a D-dimensional infinite space. Their kinetic 
energy is given by Ylii=o N-iPf / ^'^'^ they interact with each other 
via a two-body Coulombian interaction eiej/fij. We will assume that the 
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masses rrii and the charges allow the existence of a bound state. Once 
the free motion of the center of mass is discarded, we are led to a D{N — 1)- 
dimensional configuration space that can be described by the relative po- 
sitions q = {ro,j}j=i...7V-i with respect to one distinguished particle. The 
Hamiltonian is given by 

N-l N-1 N-l 

1=1 ' *j=i *j=o ■' 

The notation mjj stands for the reduced mass miirij / {nii + nij) . For D ^ 2, 
one can eliminate the Coulombian simple poles rjj = ||ro,i — rojH = in 
the local energy by choosing the trial function as follows: 

ip{q) = exp ( - - ^ Kj'^ij) (9) 
i,j=0 

with Xij = —2mijeiej/{D — 1). When this choice does not provide a nor- 
malisable function, it should be understood that the exponent is just the 
first order of a Taylor expansion near rij = 0. Whenever Q is actually 
square integrable on Q, the local energy reads 

^--E^-E^^oXi). (10) 

i<] 

The last sum involves all the angles i. A; that can be formed with all the 
triangles made of three distincts particles. This expression treats all the 
particles on an equal footing and it is clear that the local energy is bounded 
everywhere. Bounding simultaneously all the cos{j,i,k) by ±1 can be a 
rather crude approximation. One should instead take into account the cor- 
relation between the angles. For instance, when D = 3, as soon as ^ 3, 
their number (N{N — 1){N — 2)/2) exceeds the number of independent vari- 
ables minus three Euler angles and a dilatation factor 3(A'^ — 1) — 3 — 1). 
For N = 2, we recover the exact ground state since the last sum is absent 
and the local energy is constant. For a helium like atom (D = 3, = 3) of 
charge {Z — 2)e with the nucleus considered as infinitely massive compared 
to the electron mass, only two angular terms survive and we get (in atomic 
units) ii^ioc = —Z'^ — 1/4 + Z{cos9i + cos02)/2. The two angles are taken 
at the vertices made by the two electrons: the sum of their cosines is bound 
from below by and from above by 2 (diametrically opposed electrons). 
The lower bound of the local energy is -Z^- 1/4; this is not an interesting 
piece of information since we know that Eq is larger than the energy —Z'^ 
obtained by neglecting the strictly positive repulsion of the electrons. On 
the other hand the bound Eq ^ —{Z — 1/2)'^ provides a simple, analytical, 
non-trivial result. 
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5 The hydrogen atom in a magnetic field 

The third appHcation of the differential method to a non-integrable system 
will concern the hydrogen atom in a constant and uniform magnetic field 
Buz- While the positivity of the ground state wavefunction was guaran- 
teed by the so-called Krein-Rutman theorem in all the previous examples 
(applicable for any Hamiltonian of the form Q, see for instance (Reed 

6 Simon, 1978, XIII. 12)), it is not valid any longer for an arbitrary po- 
tential when a magnetic field is present (Helffer, Hoffmann-Ostenhof 
& Owen, 1999). Nevertheless, for the hydrogen atom, the attractive in- 
teraction between the nucleus and the electron keeps the orbital momen- 
tum Lz of the ground state at zero for any arbitrary value of B (Avron, 
Herbst &: Simon, 1977) and the Krein-Rutman theorem applies when re- 
stricted to the Lz = subspace. In order to preserve the symmetry of 
the ground state, the trial function will be chosen to be strictly positive 
and even with respect to z; hence we can work in the half space where 
z ^ 0. With a vanishing paramagnetic term, (p depends on the coordinates 
q = {p,r) only (see FIG. In atomic units, the local energy is given by 
-^loc ~ V — Aip/{2ip) where the effective potential is V{p,r) = B'^p^ /8 — l/r. 
In order to eliminate the Coulomb singularity, one must impose some lo- 
cal conditions on the logarithmic derivatives of ip. More precisely, with 
S = In 99, we must have drS{0,0) = — 1 and dpS{0,r) = for all r ^ 0. 
Assuming that S is smooth enough near p = 0, it takes the general form 
S{p,r) = —r + r'^l{r) + p'^h{p,r) where I and h are two smooth functions. 
Choosing / = and h = —B/A [resp. h = Q] will bound from above [resp. 
below] the local energy: —1/2 ^ Eq ^ ~l/2 + B/2. Like in the previous 
example the lower bound is useless since it can be guessed from the very 
beginning. The real challenge here is to improve the lower bound without 
introducing a divergence as r — > 00 in any direction characterized by a. Af- 
ter a detailed examination of the possible balance between the asymptotic 
behaviour of / and /i as r — > cxd, a trial function can be constructed in order 
to improve the trivial lower bound for B large enough. Namely if we take 

5 = —r — Bp'^/A + p'^{r — \J r'^ — p^) /{p^ + br/^/B) we improve the trivial 
lower bound for B > 2.3 (see FIG^. Note that for very large B we recover 
a wavefunction that mimics a Landau state. 

6 A local algorithm for improving the bounds ; 
application to the quartic oscillator 

Assume now that c/jq = exp(S'o) bounds the local energy. Is there any sys- 
tematic strategy to improve the bounds and, one day, compete with the 
very high precision of the secular variational and perturbative methods or 
numerical diagonalisation of truncated matrices ? Of course one can always 
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combine these four approaches but let us look, for the moment, how the local 
character of the differential method can be exploited further. Suppose that 
qq is a point where q ^ -E'loc'l?) reaches its lowest non-degenerate value. 
Among all the possible infinitesimal perturbations of Sq, only those that are 
localised in the neighbourhood of qo are relevant since adding a perturba- 
tion far away from qo will not affect the absolute minimum. The appropriate 
framework for local studies in an infinite functional space is bifurcation the- 
ory. Since the local energy, the determination of the critical points and their 
stability involve a finite number of derivatives, we expect that the number of 
relevant control parameters A remains finite for low-dimension configuration 
space very much like the central result of catastrophe theory (POSTON & 
Stewart, 1978; Dem azure, 2000). We wih leave this quantitative study 
for future investigations. For the moment, let us keep the discussion at a 
qualitative level only with a Id Hamiltonian of the form H = —A/2 + V and 
take S = Sq + SS with a Gaussian perturbation 5S{q) = s ex];){—{q — a)'^ /a'^) 
controlled by three parameters A = (s, a, a). The specific choice of the form 
of SS is not important here ; only a finite number of pointwise derivatives 
will matter for locally improving the bounds as long as SS does not change 
the normalisability of the trial function. The choice of a gaussian is partic- 
ularly simple: it will modify the local energy in a neighbourhood of a whose 
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size is governed by a and by the magnitude s. This perturbation is quah- 
tatively reproduced in FIGEJ The value of £"100 at g = a is increased (resp. 
decreased) for a small but finite positive (resp. negative) s. We can apply 
this procedure near the absolute minimum (resp. maximum) of Ei^c and 
repeat it for the possible absolute extrema that may have emerged during 
the previous step. We get an iteration sequence that may systematically 
improve the bounds. Still, this algorithm is slowed down because if we try 
to "lift up the dress" too much, a "prudish censor" lowers it on the both 
sides oi q = a (this phenomenon is not specific to Id). 




Figure 2: Increasing (dashed line) with one local perturbation is lim- 

ited. When adding to ln((/Jo), say, a gaussian perturbation localised near 
an isolated minimum of the local energy, we can increase the value of the 
absolute minimum by a finite amount (solid line in the left figure). But if 
the magnitude of the Gaussian is increased too much, the minimums created 
by passing through a bifurcation can decrease below the original minimum 
(solid line in the right figure). 

To be more precise, consider a quartic potential given by V{q) = r'^q'^{q^+ 
77 <5^)/2 where 77 = ±1. In order to bound the local energy as |(?| 00, we can 
use a JWKB-Iike expansion for Sq. If we want a uniformly smooth expression, 
we can take 

S,{q) = -\r{q^ + 5^'^ + \r5\l - r^){q' + 6^)'/' 

The last term is chosen to improve the trivial lower bound given by the 
minimum of V. For the arbitrary choice r = l/\/2, rj = —1 and 6"^ = 8, 
FIG.Elshows how the lower bound can be improved from —3.27 up to —2.74 
(the exact result is -2.66) when adding to Sq enough Gaussians that are equi- 
spaced by .5 with a fixed a = 1. Only their magnitude s are numerically 
optimized here, one after another. One can see a second advantage of the 
differential method when numerically implemented: not only no integral is 
required but also, provided we keep under control the instabilities that are 
illustrated in FIG. El the optimization algorithm concerns a small number 
of parameters at each step (just one in the example given in FIG. ^ to be 
compared with the large number of parameters to be optimized at one go in 
the final step of the variationnal method. 
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Figure 3: Adding n Gaussians to Sq allows to increase the minimum of 

(n) 

the local energy when constructing the sequence E^^^^. For each iteration, 
only one scalar parameter is optimized (the amplitude of the gaussian being 
added) . 

7 Conclusion 

The differential method appears to be new kind of general theoretical tool for 
obtaining rigourous information on a ground state energy. Its local character 
makes it quite different from the traditional ones (to put it succintly, the 
variational, the perturbation and the numerical diagonalization techniques). 
In this paper, I have given some qualitative and quantitative arguments to 
show how simple and efficient it can be. However, I should insist that even 
in the cases where the variational or perturbative techniques can be applied, 
the aim of the paper is not to seek for performance: for the moment the 
differential method is too young to compete by itself with the traditionnal 
methods. One short term possibility is to calculate the extrema of the local 
energy constructed with the trial function given by the other methods. The 
idea of locally modifying the local energy or any local function of the same 
type — for instance, those currently used in Monte-Carlo methods — may be 
fruitful as well (Caffarel, 2004). Actually, in order to convince the reader 
that the method is indeed applicable in a wide field of physics and furnishes 
some reasonable results, I had to compare them with some more precise ones 
and therefore I dealt with situations where the exact ground state energy 
was already known with the help of other methods. The algorithm that 
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is presented in section El is chosen to prove how the local sensitivity of the 
local energy can be exploited to systematically improve the bounds. The 
feasibility is in itself not obvious and is worth to be demonstrated even in 
the simplest cases. 

This work could not have been started without Hector Giacomini's bril- 
liant intuition that some relevant information on Eq could be extracted 
from ((21). I am very indebted to Dominique Delande and Benoit Gremaud 
for sharing their penetrating thoughts, their skilled numerical calculations 
in Coulombian problems and, not least, their kindful hospitality at the Lab- 
oratoire Kastler Brossel. 



References 

Arnold V.I. (1984): Catastrophe Theory. New York: Springer- Verlag. 

AvRON J., Herbst I. & Simon B. (1977): "The Zeeman effect revisited", 
Phys. Lett. A, 62, pp. 214-216. 

Barnsley M. (1978): "Lower bounds for quantum mechanical energy lev- 
els", J. Phys. ^, 11(1), pp. 55-68. 

Barta J. (1937): "Sur la vibration fondamentale d'une membrane", C. R. 
Acad. Sci. Paris, 204(7), pp. 472-473, (in french). 

Baumgartner B. (1979): "A class of lower bounds for Hamiltonian oper- 
ators", J. Phys. A, 12(4), pp. 459-467. 

Bessa G.P. & Montenegro J.F. (2004): "An Extension of Barta's The- 
orem and Geometric Applications", eprint arXiv:math/0308099 

BoHiGAS O., Boose D., Egydio de Carvalho R. & Marvulle V. 
(1993): "Quantum tunneling and chaotic dynamics". Nuclear Phys. A, 
560, pp. 197-210. 

Caffarel M. (2004): private communication. 

Crandall R.E. & Reno M.H. (1982): "Ground state energy bounds for 
potentials ||x||''", J. Math. Phys., 23(1), pp. 64-70, erratum[J. Math. 
Phys. 23, 1737 (1982)]. 

Demazure M. (2000): Bifurcations and Catastrophes, Universitext. 
Springer. 

DuFFiN R.J. (1947): "Lower bounds for eigenvalues", Phys. Rev., 71, pp. 
827-828. 



12 



Helffer B., Hoffmann-Ostenhof T. & Owen M. (1999): "Nodal sets 
for groundstates of the SchrOdingcr with zero magnetic fields in a non 
simply connected domain", Comm. Math. Phys., 202, pp. 629-649. 

Maslov v. p. & Fedoriuk M.V. (1981): Semi- Classical Approximation 
in Quantum Mechanics, vol. 7 of Mathematical physics and applied math- 
ematics. Dordrecht: D. Reidel publishing company. 

POSTON T. & Stewart I. (1978): Catastrophe Theory and its Applications. 
London: Pitman. 

Reed M. & Simon B. (1978): Analysis of operators, vol. 4 of Methods of 
modern mathematical physics. New York: Academic Press. 

SCHMUTZ M. (1985): "The factorization method and ground state energy 
bounds", Phys. Lett. A, 108(4), pp. 195-196. 

Thirring W. (1979): Quantum mechanics of atoms and m,olecules, vol. 3 
of A course in mathematical physics. New York: Springer- Verlag. 



13 



